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Coupling Index and stocks 
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Abstract 



o 
o 

£N) ■ In this paper, we are interested in continuous time models in which the index level induces some 

feedback on the dynamics of its composing stocks. More precisely, we propose a model in which the 
log-returns of each stock may be decomposed into a systemic part proportional to the log-returns of the 
y-^ • index plus an idiosyncratic part. We show that, when the number of stocks in the index is large, this 

model may be approximated by a local volatility model for the index and a stochastic volatility model 
, for each stock with volatility driven by the index. This result is useful in a calibration perspective : it 

suggests that one should first calibrate the local volatility of the index and then calibrate the dynamics 
of each stock. We explain how to do so in the limiting simplified model and in the original model. 
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Introduction 

> 

■ From the early eighties, when trading on stock index was introduced, quantitative finance faced the 

^) \ problem of efficiently pricing and hedging index options along with their underlying components. Many 

advances have been made for single stock modeling and a variety of solutions to escape from the very 
restrictive Black & Scholes model has been deeply investigated (such as local volatility models, models with 
jumps or stochastic volatility models). However, when the number of underlyings is large, index option 
pricing, or more generally basket option pricing, remains a challenge unless one simply assumes constantly 
correlated dynamics for the stocks. The problem then is the impossibility of fitting both the stocks and the 
index smiles. 

We try to address this issue by making the dynamics of the stocks depend on the index. The natural 
fact that the volatility of the index is related to the volatilities of its underlying components has already 
been accounted for in the works of Avellaneda et al. pQ and Lee et al. [13]. In the first paper, the authors 
use a large deviation asymptotics valid for small values of the product of the maturity by the square of the 
volatility to reconstruct the local volatility of the index from the local volatilities of the stocks. They express 
this dependence in terms of the implied volatilities using the results of Berestecky et a/. ([I] In the second 
paper, the authors reconstruct the Gram-Charlier expansion of the probability density of the index from the 
stocks using a moments- matching technique. Both papers consider local volatility models for the stocks and 
a constant correlation matrix but the generalization to stochastic volatility models or to varying correlation 
coefficients is not straightforward. 

Another point of view is to say that the volatility of a composing stock should be related to the index 
level, or say to the volatility of the index, in some way. This is not astonishing since the index represents 
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the move of the market and reflects the view of the investors on the state of the economy. Moreover, it is 
consistent with equilibrium economic models like CAPM. Following this idea, we propose a new modeling 
framework in which the volatility of the index and the volatilities of the stocks are related. We show that, 
when the number of underlying stocks tends to infinity, our model reduces to a local volatility model for the 
index and to a stochastic volatility model for the stocks where the stochastic volatility depends on the index 
level. This asymptotics is reasonable since the number of stocks composing an index is usually large. As a 
consequence, the correlation matrix between the stocks in our model is not constant but stochastic and we 
show that it is consistent with empirical studies. Finally, we address calibration issues and we show that it 
is possible, within our framework, to fit both index and stocks smiles. The method we introduce is based 
on the simulation of SDEs nonlinear in the sense of McKean, and non-parametric estimation of conditional 
expectations. 

This paper is organized as follows. In Section 1, we specify our model for the index and its composing 
stocks and in Section 2 we study the limiting model when the number of underlying stocks goes to infinity. 
Section 3 is devoted to calibration issues. Numerical results are presented in Section 4 and the conclusion is 
given in Section 5. 

Acknowledgements: We thank Lorenzo Bergomi, Julien Guyon and all the equity quantitative research team 
of Societe Generale CIB for numerous fruitful discussions and for providing us with the market data. 



1 Model Specification 

An index is a collection of stocks that reflects the performance of a whole stock market or a specific sector 
of a market. It is valued as a weighted sum of the value of its underlying components. More precisely, if if 1 
stands for the value at time t of an index composed of M underlyings, then 

M 

/f = 5>^ M , (1) 

i=i 

where S{' M is the value of the stock j at time t and the weightings (ti>j)j=i...M ar e given constant^). 

Unless otherwise stated, we always work under a risk-neutral probability measure. In order to account 
for the influence of the index on its underlying components, we specify the following stochastic differential 
equations for the stocks 

Vj G {1, . . . , M}, -j— = (r - 6j)dt + a(t, I t M )dB t + Vj (t, S>' M )dWl , Sq = 4 (2) 

where 

• r is the short interest rate, 

• Sq is the initial value of the stock j, 

• 8 3 ■ G [0, oo[ is the continuous dividend rate of the stock j, 

• (3j is the usual beta coefficient of the stock j that quantifies the sensitivity of the stock returns to the 
index returns (see the seminal paper of Sharpe [17j). It is defined as ^y^'^ where (respectively 
rj) is the rate of return of the stock j (respectively of the index), 

2 The weightings are periodically updated but, as usually assumed, we suppose that, up to maturities of the options considered, 
they are constant. When updated, they are often chosen proportional to the market capitalizations of the stocks. 
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(Bt)te[o,T], {Wt)te[o,T],- ■ ■ , (W* M )te[0,T] are independent Brownian motions, 



• the functions a, r)i, . . . ,t]m : [0,T] x R — > R satisfy the usual Lipschitz and growth assumptions that 
ensure existence and strong uniqueness of the solutions (see for example Theorem 5.2.9 of [12] ) : 



(HI) 3K such that V(i, s u s 2 ) G [0, T] x R x R 



M 
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M 

E 
E 
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\ fc=l / 
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V fe=i / 
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fe=l 



< K\ Sl -s 2 \ 



< K\si-s 2 \ 



As a consequence, the index satisfies the following stochastic differential equation 



M 

dI t M = rI t M dt-[Y / S j w :i S l 
^=1 




*{t,I?)dB: 



M 



z ^w J Si' M Vj (t,S(- M )dW^ 

3=1 



Before going any further, let us make some preliminary remarks on this framework. 



(3) 



We have M coupled stochastic differential equations. The dynamics of a given stock depends on all 
the other stocks composing the index through the volatility term a(t, / t M ). Since there are M linearly 
independent assets and M + 1 driving Brownian motions, the market is incomplete. 

Accounting for the dividends is not relevant for all types of indices. Indeed, for many performance- 
based indices (such as the German DAX index) dividends and other events are rolled into the final 
value of the index. 



The cross-correlations between stocks are not constant but stochastic 

pij (t) 



Note that they depend not only on the stocks but also on the index. More importantly, it is commonly 
observed that the more the market is volatile, the more the stocks tend to be highly correlated. This 
feature is reproduced here as we can easily check that an increase in the index volatility, with everything 
else left unchanged, produces an increase in the cross-correlations. 

In a recent paper, Cizeau et al. [8] show that it is possible to capture the essential features of stocks cross- 
correlations, in particular in extreme market conditions, by a simple non-Gaussian one factor model. 
The authors successfully compare different empirical measures of correlation with the prediction of the 
following model : 

r j (t)=/3 j r I (t)+e j (t) (4) 



where rj(t) = —f 1 is the daily return of stock j, rj(t) is the daily return of the market and the 

n 

residuals tj (t) are independent random variables following a fat-tailed distributioro 



j The authors have chosen a Student distribution in their numerical experiments. 
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Our model is in line with J3]). Indeed, since the beta coefficients are usually narrowly distributed 
around 1, the factor J2j=i ftj w jSt' M °f < T (^^t M ) m © 1S close to I t M . Moreover, since 

(/ r T M \ 2 \ M „ T M 

I / J2 w i S i' M vAt,S? M )dWn U^^jf E S*"))')* ~X>^, 

one can neglect the term Sjli w j^t' M Vj(t> Sf' M )dW^ in the dynamics of the index when Ylj=i w j is 
small. Of course, this approximation worsens when the maturity T increases. The latter condition is 
satisfied when M is large and the weighting vector (wi, . . . , wm) is close to the vector (A, . . . , A) with 

constant coefficients for which Ylj=i IT 7 = ~M' Then, if we denote by rj the log-return of the stock j 
and by fjM the log-return of the index, both on a daily basis, we will have 

r,j = (3jTiM + VjAW 3 + drift, 

where is an independent Gaussian noise. Consequently, in our model too, the return of a stock 

is decomposed into a systemic part driven by the index, which represents the market, and a residual 
part. 



2 Asymptotics for a large number of underlying stocks 

The number of underlying components of an index is usually largfl As discussed in the previous section, 
when YljLi w j is small, one can neglect the term Ylj=i w jSt' Mr lj(t> S{ ,M )dW^ in ([3]) and derive a simplified 
approximate dynamics for the index. The aim of this section is to quantify the error we commit by doing so. 

To be specific, consider the limit candidate (It)te[o,T\ solution of the following SDE : 

f dl t = (r-5)I t dt + [3I t a{t,I t )dB t 

I Jo = io (5) 

where io = Sj=i w j s o an< ^ ^ an d P are t w0 constant parameters that will be discussed later. 

In the following theorem, we give an upper bound for the L 2p -distance between (^t M )tg[o,T] an d (Jt)te[o,ri 
under mild assumption on the volatility coefficients : 



Theorem 1 — Let p E N*. Under assumption fTlQP and if the following assumptions on the volatility 
coefficients hold, 

(H2) 3K b such that V(t,s) € [0,T] X R + , \a(t, s)\ + \r)j(t, s)\ < K b . 

(H3) 3K a such thatV(t,si,s 2 ) € [0,T] x R + x R + , \ Sl a(t, Sl ) - s 2 a(t,s 2 )\ < K«\s x - s 2 \. 
then 



E 



( sup \i? -i t A<c T \ ff; 



M 



2p / \ 2p\ 

' M - ' 



where 



and 



C T = ^-^{TP + K p Kl P )C p exp (4 2p - 1 T(2 2p - 1 K p T p - 1 (f3K a ) 2p + (2T) 2p - 1 S 2p + r 2p T 2p - 1 )) 
C p = ^mx^ \s 3 \ 2p exp ^ \2r + (2p - l)(max^ + l)K 2 ^J pT^j . 



4 500 stocks for the S&P 500 index, 100 stocks for the FTSE 100 index, 40 stocks for the CAC40 index, etc. 
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According to this result proved in the appendix, the smaller P^f = 

... I r n i.i i t 1 T M 



pM 

r 



and Pg 1 = f X)j=i ^j'l^j ~~ the closer / and I M . The first quantity P,^ is small when the weighting vector 
(u>i, . . . , wm) is close to (i, . . . , and M is large. Let us now discuss how to choose /? and S minimizing 



Pf and P, 



Then 



A/ 



Let Yg and Y$ be discrete random variables having the following probability distributions 
Vj G {1, . . . , M}, P (Yp = pj) = J' and P (Y$ = <>, 



EM 
1=1 



pM 




M 



xE\Y f3 -f3\ and 



pM 
^8 



Wi x E \Y 5 - S\ 



Consequently, the optimal choice of the parameters is the median^ of Yp for p and the median of Yg for 

for the simplified 

index dynamics, one should take (3 



piv-uciuiuii u± a,o Var(ri) ""^^ io v^LLcm. uu 

1. In Table [l] we sec that on the example of the Eurostoxx index at 



December 21 2007, the optimal choice of (3 is very close to 1 and that the quantities of interest, (Pp^ pt ) 2 and 
(PfiLi) 2 are also very close to each other. 





Popt 






0.026 


0.975 


0.0173 


0.0174 



Table 1: Computation of (P™) 2 , / 9 opt and {Pp ) for the Eurostoxx index at December 21, 2007. The beta 
coefficients are estimated on a two year history. 

The next theorem states that, under an additional assumption on the volatility coefficients, the L 2p - 
distance between a stock (£'t ' M )te[o,T] and the solution of the SDE obtained by replacing I M by I 



Si 



(r - 6j)dt + /3j a(t, I t )dB t + ^(t, S 3 t )dWl 



I 1 
d 



(6) 



is also controlled by 2p-powers of P^f , P¥ and Pf 1 ■ One major drawback of the limiting simplified model ([5])- 
([6]) is that the limit index It is only approximately equal to the reconstructed index level I t d = Ejf=i w j^t ■ 
The next result also gives an estimation of the difference between I M and I M in terms of P* i , P^' 1 and P/ f , 

—M 

which combined with the previous theorem, provides an estimation of the difference between I and / . 



Theorem 2 — Let p e N*. Under the assumptions of Theorem 1 and if 

(HA) 3K V such thatVj, V(t,Si,s 2 ) G [0,T] x R + x R + , [siTf^t, s x ) - s 2 T]j(t, s 2 )\ <K n \s 1 -s 2 \- 

3K Lip such that V(i, Si,s 2 ) G [0,T] x P + x R + , \a(t,si) - a(t, s 2 )\ < K Lip \si - s 2 \. 
Then, Vj G {1,...,M} 7 

2p / \ 2p\ 



E sup \Si' M -S{\ 2 *> 

\0<t<T 



5 The median of a real random variable X is any real number m satisfying : 

P(X < m) > i and P(X > m) > -. 

It has the property of minimizing the L 1 -distance to X : m = argminElX — x\. 

xeK 
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where 

C 3 T = Q 2p ~ 1 K p T P (3 2p C% Kl P e 3 2p ' 1 ({r-S J fT 2 "- 1 +K p T'- 1 K 2 rl ''+2^- 1 K p T'- 1 l3^K 2 b p )T^ 



3 Lip 

Moreover, for I t = Y^,jLi w jSt> one has 

2p\ 



where Ct = max CL. 

l<j<M 



The proof can also be found in the appendix. In the following corollary, we consider the limit M — > oo 
supposing that the weight of the j-th stock, now denoted by Wj 1 , depends on M. 



Corollary 3 — Under the assumptions of Theorems 1 and 2 and if 

(H5) there exists a finite constant A such that max I (sq) 2 + (ftj) 2 + (Sj) 2 ) < A, 

j>i \ J 



M 

(H6) 7 M = J2 <4 -r^ *o e (0, +oo) ; 



(H7) P™ = 



M 

. V« f ) 2 — > o, 

\ 3 = 1 



M 

(H8) P M = J2wf\P j -[3\-^ 0, 



M 



t/ien, /or any p G N*, one /las 

0. 



E^sup |/ t M -7 t | 2p> ) — ► andVjeN*, E f sup |5f M -5|'| 2p 

\0<t<T J M ~> 00 \0<t<T 

//, m addition, sup^^ wf 1 < oo, i/iera E ^sup 0<t<T |/ t M — /f f | 2p 



m . , 



Assumptions (7lH]), (71©, (7"4H]) an d hold for instance when u>* f = for 1 < j < M and — >j^oo 
/9j — >j^oo P and <5j — ^ 
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Simplified model 

To sum up, we have shown that, under mild assumptions, when the number of underlying stocks is large, 
the original model may be approximated by the following dynamics 

dS j 

Vj G {1, . . . , M}, —f = (r- Sj)dt + 13, a(t, I t )dB t + Vj (t, S 3 t )dWi 

(7) 

^ - {r-8i)dt + a{t,I t )dB t . 
n 

Of course the distance between this limiting model and the original one increases with the maturity. 

Interestingly, we end up with a local volatility model for the index and, for each stock, a stochastic 
volatility model decomposed into a systemic part driven by the index level and an intrinsic part. The 
calibration procedures presented in the next section are based on this intuition : even in the original model, 
we are going to calibrate a as if it was the local volatility function of the index. 

Note that this simplified model is not valid for options written on the index together with all its composing 
stocks since the index is no longer an exact, but an approximate, weighted sum of the stocks. In this case, 
one should consider the reconstructed index I t — X^fci ^j^i or use the original model. The simplified 
model can be used for options written on the stocks or on the index or even on the index together with few 
stocks. 



3 Model calibration 

Calibration, which is how to determine the model parameters in order to fit market prices at best, is of 
paramount importance in practice. In the following, we try to tackle this issue for both our simplified and 
original models. 

3.1 Simplified model 

In the simplified limiting model, the only factor which influences the dynamics of a given stock is the 
simplified index I t which evolves according to an autonomous SDE. So it is enough to address the calibration 
of a given stock together with the index and we drop the index j of the stock for notational simplicity. 

S)dt + (3a(t, I t )dB t + n(t, S t )dW t , S = s 

(8) 

5i)dt + a(t,I t )dB t , I = i . 

The short interest rate and the dividend yields can be extracted from the market. The calibration of the local 
volatility function a to fit index option prices is a classical problem. According to Dupire [10], if Ci(t,K) 
denotes the market price of the call option with maturity t and strike K written on the index, then for 

2, R) = 2 j£& K ) + ( r - WfjK*. K) + SjCijt, K) 

one has Cj(t, K) — E (e~ rt (It — K) + ) for all t, K > 0. Of course, in practice the market quotes call options 
only for a finite number of couples (t, K). What seems a common practice among banks is to look for a in a 
parametric family of functions and compute the parameters minimizing the distance between these quoted 
prices and the call prices associated with the parametrized local volatility function. Since each practitioner 
may choose his favorite procedure to address this classical problem of local volatility calibration, we will not 
enter in more details. 



dS t 



= (r- 



dl t , 
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We also assume that a local volatility function is associated with the stock by the same procedure and 
denote by vi oc (t,x) the local variance function of the stock computed as the square of this local volatility 
function. So the local volatility model 

^l = ( r -6)dt+ ^v loc {t,S t )dWt ,So = s Q (9) 
St 

is calibrated to the quoted prices of vanilla options written on the stock. In |8"]). by independence between 
B and W, the variance of the stock at time t is equal to P 2 a 2 (t, I t ) + rf(t, St). According to Gyongy [TT], if 

Vt,x>0, E(f3 2 a 2 {t,I t )+T 1 2 {t,S t )\S t = x) =v loc (t,x) 

then ijHJ) and the local volatility model ((9]) induce the same marginal distributions for the stock and therefore 
the same prices for the vanilla call options written on it : IE (e~ rt (St — K) + ) = E (e~ rt (St — K) + ) for all 
t, K > 0. Hence if 

Vi, x > 0, r)(t, x) = ^v loc (t,x)~ p 2 £{a 2 (t,I t )\S t = x), (10) 

then the stock dynamics in ||5J) is calibrated to the quoted prices of the vanilla options written on the stock. It 
remains to choose the coefficient (3 and the function -q so that this equality is satisfied. The fact that the law 
of (St, It) given by © and therefore the conditional expectation in (|10[) depend on (P, -q) makes this problem 
difficult. Nevertheless, intuitively, when one fixes a value of p that is not too large, one should be able to find 
a function r\ such that (fTtl)) is satisfied. The calibration of the stock smile seems over-parametrized and one 
should rely on the interpretation of /3 as a regression coefficient to choose its value. This issue is discussed 
in the next section. Then we explain how to approximate the conditional expectation and deduce r\ for a 
fixed value of (3. 

Let us already point out that the calibration of our simplified model gives an advantage to the fit of index 
option prices in comparison with options written on the stocks, which is in line with the market since index 
options are usually very liquid in comparison with individual stock options. 

3.1.1 Choice of the coefficient (3 

The interpretation of j3 as the regression coefficient of the log-returns of the stock with respect to the 
log-returns of the index makes it possible to estimate this coefficient on historical data. Nevertheless, when 
the historical estimator fihist is large, then the difference in the r.h.s. of (TlU|) may become negative for some 
(t, x) when p = Phist ■ Then the square root is no longer defined and calibration for this choice of (3 is no 
longer possible. 

In Figure [TJ we have plotted the local volatility of the stock x i— > ^vi oc (T, sqx), the local volatility 
of the index x i— > a(T,iox), the systemic part of the volatility of the stock x i— * (3histc(T, sqx) and x i— > 
Phiat a/E (<t 2 (T, It)\St = sqx) when r\ is set to zero (which intuitively gives the lowest local volatility function 
of the stock that one can obtain in our model (O) as functions of the moneyness for a maturity T = 1 year. 
We considered three representative components of the Eurostoxx which is composed of M = 50 stocks : 
AXA, ALCATEL and CARREFOUR at December 21, 2007. We made this choice deliberately in order to 
point out the extreme situations that one can face : 

• AXA is an example of a stock with a high historical beta coefficient (Phist = 1-4), 

• CARREFOUR is an example of a stock with a low historical beta coefficient (Ptist = 0.7), 

• ALCATEL is an example of a stock with a high volatility level but with a rather flat smile (Phist = 1-1). 

Clearly, we can deduce that the market is choosing a P coefficient for both AXA and ALCATEL that 
is lower than the historical one whereas, for CARREFOUR, one can plug the historical P, or even a larger 
one, in ([5]) and still be able to calibrate the model. 
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AXA 



9 <7(T, i x) 
"+ /3ct(T, tQi) 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.1 

Moneyness 



® <7(T, i x) 
"+ /3ct(T, i i) 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.1 



CARREFOUR 




_^ ; --A- 



0.0 0.2 0.4 



0.8 1.0 1.2 

Moneyness 



Figure 1: Local volatility x i— > y/ vi oc (T, sqx) together with x i— > a(T,iox), x i— » PhistO'(T,iox) and cc i— > 
/3/ MS tv / E(cr 2 (T,i T )|S , T = s x) when rj is set to zero for AXA, ALCATEL and CARREFOUR. 



A satisfactory way to handle the estimation of the beta coefficient would be to compute an implied beta 
calibrated to the prices of options involving the correlation between the stock and the index. Unfortunately, 
no such option is liquid in the market (the most liquid correlation swaps are sensitive to an average correlation 
between all the stocks composing the index). 

So we suggest to choose 

a ■ in ■ r V v ioc(t, s x)\ 

(3 = mm PusU mf - — TT- — % — ■ (H) 
\ t,x>o a(t,ioX) I 

Even if we have no proof that this choice of beta makes the calibration possible, it is sensible and we have 
checked that it works on the three examples of AXA, ALCATEL and CARREFOUR. 

When one is interested in options written on the index together with all its components, one should use 

the reconstructed index level I t = Ylj=i w j^t instead of I t . Of course, the reconstructed index dynamics 

will reproduce the quoted prices of vanilla options written on the index all the better as is close to the 
calibrated limiting index level It- According to Theorems Q] and [2l for this latter property to hold, one 
needs PpL 1 — J2j=i w i\Pj ~ M — I Hj=i w jPj ~ w j\ to ^e small. When, because of the minimum in 

equality (fTT|) . PaLi is larger for the actual choice of coefficients (3 than for the historical choice, one may 
take larger values of beta for stocks like CARREFOUR to decrease PgL-y and improve the calibration of the 
reconstructed index. 

3.1.2 Estimation of the conditional expectation 

The idea behind the following techniques is to circumvent the difficulty of calibrating the volatility 
coefficient r\. Indeed, if we plug the formula (| 10[) in ([8]), we obtain a stochastic differential equation that is 
nonlinear in the sense of McKean : 

J Q 

-± = (r - 5)dt + (3 a(t, I t )dB t + ^v loc {t, St) - (3 2 E(a 2 (t, I t ) \ S t )dW t , S = s 

dl (12) 
_* = (r - 6 T )dt + a{t,I t )dB u I = i . 
h 

For an introduction to the stochastic differential equations nonlinear in the sense of McKean and to 
propagation of chaos, we refer to the lecture notes of Sznitman [18] and Meleard |15) . In our case, the 
nonlincarity appears in the diffusion coefficient through the conditional expectation term. This makes 
the natural question of existence and uniqueness of a solution very difficult to handle. The case of a 
drift coefficient involving a conditional expectation has only been handled recently even for a constant 
diffusion coefficient (see for instance Talay and Vaillant [12] and Dermoune 0). Meanwhile, it is possible 
to simulate such a stochastic differential equation by means of a system of N interacting paths using either 
a non-parametric estimation of the conditional expectation or regression techniques. The advantage of the 
regression approach over the non-parametric estimation is that it also yields a smooth approximation of 
the function x i— > E (<r 2 (t, It ) \ St = x) whereas, with a non-parametric method, one has to interpolate the 
estimated function and to carefully tune the window parameter to obtain a smooth approximation. 

3.1.2a Non-parametric estimation 

Non-parametric estimators of the conditional expectation, and more generally non-parametric density 
estimators, have been widely studied in the literature. We will focus on kernel estimators of the Nadaraya- 
Watson type (see [3D] and [TDJ) : given N observations {S^t, Ii,t)i=i...N of (S t ,It), we consider the kernel 
conditional expectation estimator of E (er 2 (t, It) \ S t = x) given by 
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where if is a non- negative kernel such that J R K(x)dx = 1 and Tijv is a smoothing parameter which tends 
to zero as N — > +00. This leads to the following system with N interacting particles : VI < i < N, 



dS it 



(r - 5)dt + /3a(t,I iit )dB iit + J v loc (t, S itN , t ) - P 



,N,t-^j,N,t 



-dWij, S itNt0 = s 



-jM. = (r - + I iA )dB iA , 7 i>0 = i 



(13) 

where (JS^, Wi)i>i is a sequence of independent two-dimensional Brownian motions. The integer i indexes 
the sample-paths of the fixed stock that we consider. In their dynamics, the conditional expectation term 
has been replaced by interaction. The price in the calibrated model of a European option with maturity T 
and payoff function h : C([0, T], R) — > R written on the stock may be approximated by 



1 N 

-J2e- rT h(S i;N ,.). 



(14) 



The 2A-dimensional SDE may be discretized using the Euler scheme. Let n S N* and = to < • ■ ■ < 
t n = T be the subdivision with step £ of [0, T]. For each k G {0, . . . , n - 1}, V 1 < i < N, 



Si 



N,t k 



Si 



\ 



Vloc(tk,Si t N,t k ) - 0* 













Ef =1 k 


( S i,N 









+ (r-<5)^+/3a(t fe ,/. t Jy , ^G l , fc 



where (C?j,fe)i<i<iV,o<fc<n-i an d {Gi,k)i<i<N,o<k<n-i are independent centered and reduced Gaussian ran- 
dom variables. 



3.1.2b Parametric estimation 



Another approach to estimate conditional expectations is to use parametric estimators, or projection. 
This idea has also been widely used and studied previously (for example in finance, one can think of 
the Longstaff-Schwartz algorithm for pricing American options [14 ). Noting that the conditional ex- 
pectation is a projection operator on the space of square integrable random variables, one can approxi- 
mate E (cr 2 (t, I t ) I St = x) by the parametric estimator X);=i a ifi( x ) where (fi)i=i...L is a functional basis 
and a = (oti)i=i...L is a vector of parameters estimated by least mean squares : given N observations 

(Si,t,Ii,t)i=i...N of (S t , I t ), a minimizes J2i=i (^fo-A'.O ~ Ya=i &lfl(Si,t) ) ■ 



3.1.3 Numerical results 
3.1.3a A toy example 

We try to calibrate a stock with a local variance function vi oc constant and equal to v. We choose <j as 
the local volatility function of the Eurostoxx index fitted to the market at December 21, 2007. 
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We simulate the system of N interacting paths p3[) and price call options for different strikes using (|14[) . 
In Figure [2 we plot the implied volatility at T = 1 obtained for independent simulations of N — 5000 paths 
and see that they are indeed close to the desired volatility level y/v. This example was generated with the 
following arbitrary set of parameters : 

S = 100, = 0.7, r = 0.05, 5 = Si = 0, ^ = 0.6, N = 5000, n = 20. 

In this example and for all the following numerical experiments, we use a Gaussian kernel : K{u) = 
-^=e 2 . The smoothing parameter Hm is set to N s which is the optimal bandwidth that one obtains 
when minimizing the asymptotic mean square error of the Nadaraya- Watson estimator under some regularity 
assumptions and assuming independence of the random variables involved (see for example Bosq [5]). 



Exact 
Simulated 



0.57 

0.80 



0.85 0.90 



0.95 1.00 1.05 

Moneyness 



I i I i I 
1.10 1.15 1.20 



Figure 2: Implied volatility obtained for nine independent simulations with N — 5000 paths. 



3.1.3b An example with real data 

In the following, we test our model with real data. More precisely, given the local volatilities of the 
Eurostoxx index and of Carrefour at December 21, 2007, we simulate (fT2"]) by different methods for a one 
year maturity. 

An acceleration technique 

The simulation of the particle system is very time consuming : for each discretization step and for 
each stock particle, one has to make N computations which yield a global complexity of order 0{nN 2 ) 
where n is the number of time steps in the Euler scheme. Acceleration techniques are thus desirable. One 
possible method consists in reducing the number of interactions : instead of making N computations for 
each estimation of the conditional expectation, one can neglect interactions which involve particles which are 
far away from each other. When the kernel used is non increasing with the absolute value of its argument, 
the easiest way to implement this idea is to sort the particles at each step and, whenever a contribution of 
a particle is lower than some fixed threshold, to stop the estimation of the conditional expectation. 
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Of course, by doing this, we lose in precision for the same number of interacting particles, especially for 
deep in/out of the money strikes. But what we gain in terms of computation time is much more important 
: in Figure [31 we plot the implied volatility obtained by the naive method and the method with the above 
acceleration technique for the same number N = 10000 of particles. We take as threshold -k and set 
/ijy = N~t° for the bandwidth parameter!! and n = 20 for the number of time steps in the Eulcr scheme. 
The computation time, on a computer with a 2.8 Ghz Intel Penthium 4 processor, is of 52 minutes for the 
naive method and of 5 minutes for the accelerated one. 




Figure 3: Comparison between the naive technique and the accelerated one for the same number N = 10000 
of particles. 



More importantly, we see that the implied volatility ct/v obtained by simulation of the system with TV 
interacting particles converges to the exact implied volatility o exac ± computed from quoted option prices as 
N tends to oo : see Figure [4] and Table [2] With a reasonable number of simulated paths, N = 200000, the 
error on the implied volatility remains clearly tolerable for practitioners (of the order of 10 bp) except for a 
deep in the money call (K — O.3S0) where it attains 195 bp. 

6 In order to smooth the estimation, one has to choose a bandwidth parameter that is greater than the theoretical optimal 
_ \ 

parameter N 5 . 
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Moneyness (^-) 


0.30 


0.49 


0.69 


0.79 


0.89 


0.99 


1.09 


1.19 


1.28 


1.48 


1.98 


Error : \a N - d exact \ 


195 


36 


8 


5 


2 


1 


2 


9 


17 


32 


56 



Table 2: Error (in bp) on the implied volatility with N = 200000 particles. 




Figure 4: Convergence of the implied volatility obtained with non-parametric estimation as N — > +oo. 



Independent particles 

Unlike the parametric method, non-parametric estimation of the conditional expectation gives the value 
of the intrinsic volatility rj at the simulated points only. However, using an interpolation technique, one can 
first reconstruct rj with Ni dependent particles and then simulate N2 independent paths of the 2-dimensional 
stochastic differential equation (|8|). By doing so, we speed up the simulations but one has to choose carefully 
the size Ni of the particle system in order to have a reasonable estimation of the intrinsic volatility and 
to tune the bandwidth parameter in order to smooth the estimation (our numerical tests were done with 

Ni = 1000, N 2 = 100000 and h Nl = N~ w ). In Figures E] and El we plot the local volatility function 
\J vi oc {t, x) and the intrinsic volatility function r)(t, x) of the stock. This latter is used to draw independent 
simulations of the index along with the stock and we see in Figure [7] that the implied volatility obtained is 
close to the right one, especially near the money. 
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Figure 5: Local volatility function yj vi oc (t, x) of Figure 6: Intrinsic volatility function r](t, x) of the 

the stock. stock. 




3.2 Original model 



We now turn to the calibration of our original model : 

Vj 6 {1, ... , M}, = (r- 5 3 )dt + (3, a(t, I t M )dB t + Vj (t, S^ M )dWi with if = ]T Wl S\' M . (15) 

St J= i 

It is rather complicated to have a perfect calibration for both index and stocks within this framework. 
Nevertheless, Theorem[T]ensures that the error of calibration of the index smile is small (at least when the ma- 
turity is not too large) when a is chosen as a local volatility function fitted to this smile. We also suppose that 

a local volatility function \J vj oc has been fitted to the market smile of each stock j . For the choice of the coef- 
ficients f3 j, we proceed like in Section l3.1.1l The coefficients r/j (t, x) = \Jv{ oc {t, x) — (3"jE(a 2 (t, i t M )|S'/' M = x) 
are then calibrated all at the same time using an adaptation of the non-parametric method presented above 
based on the simulation of N interacting (M + l)-dimensional paths. 

In comparison with the simplified model, we introduce in the calibration of the index a small error which 
grows with the maturity T. But we guarantee the additivity constraint J t M = X)i=i WiS\' M . Note that a 
similar error spoils the calibration of the reconstructed index in the simplified model (see the discussion at 
the end of Section I3.1.1|) . 

In what follows, we illustrate the effect of Theorems 1 and 2 and compare our models with a constant 
correlation model. 



4 Illustration of Theorems 1 and 2 and comparison with a constant 
correlation model 

The objective of this section is to compare index and individual stock smiles obtained with three different 
models : our original model (TT5")) . the simplified one (after letting M — » oo) and a model with constant 
correlation coefficient. More precisely, we consider the following dynamics 

1. The original model 

Vj e {1, • . • , M}, = rdt + a{t, lf I )dB t + ij(t, S 3 t M )dW 3 t with I t M = ^ Wl S\' M . (16) 

2. The simplified model 

dsi 



Vj G {1, . . . , M}, ^L=rdt + a(t, I t )dB t + V (t, S})dW? 
St 

-p. = rdt + a{t,I t )dB t . 
It 

Here we can also compute the reconstructed index I t = YliLi WiS l t . 
3. The "market" model 

dS 3 t 



(17) 



Vj G {1, . . . , M}, —L = rdt+ ^v loc (t, Sl)dW J t (18) 



with, Vi ^ j, d< W\ W j > t = pdt. 
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We deliberately dropped the dividend yields and the beta coefficients in order to simplify the numerical 
experiment. For the function a, we take as previously the calibrated local volatility of the Eurostoxx. For 
rj, which does not depend on j, we choose an arbitrary function of the forward moneyness and we evaluate 
vi oc such that the "market" model and the simplified model yield the same implied volatility for individual 
stocks. According to [TT], it is enough to take 

vi oc (t,x) = rf{t,x) + E(a 2 (t,I t )\Sl = x) 

where the conditional expectation is approximated using the non-parametric method presented above. 

Finally, we fix the correlation coefficient p such that the market model and the simplified one have the 
same ATM implied volatility for the index. 

The implied volatilities for the index and for an individual stock obtained by the three models are plotted 
in Figures [H and [5) We also give the difference in basis points between the implied volatilities obtained with 
the simplified model and the original one in Tables [4J [5] and [3j The parameters we use in our numerical 
experiment are the following : 

- Sq = ■ ■ ■ = S = 53, 

- M, Io and the weights W\, . . . , wm ■ the same as of the Eurostoxx index at December 21, 2007, 

- r = 0.045, 

- Maturity T = 1 year, 

- Number of time steps: n = 10, 

- Number of simulated paths : N = 100000. 
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Table 3: Difference (in bp) the implied volatilities of an individual stock obtained with the simplified model 
and with the original model. 
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Table 4: Difference (in bp) between the implied volatilities of the index obtained with the simplified model 
and with the original model. 



Moneyness (^) 


0.5 


0.8 


0.9 


0.95 


1 


1.05 


1.1 


1.2 


1.3 


1.55 


1.85 


2 


I*? reconstruct & original \ 


10 


5 


4 


3 


2 


1 


2 


5 


4 


1 









Table 5: Difference (in bp) between the implied volatility of the reconstructed index / in the simplified 
model and the implied volatility of the index in the original model. 

As suggested by Theorems 1 and 2, we see that the original model and the simplified one yield implied 
volatility curves that are very close to each other, both for the index and for individual stocks. The difference 
in basis points between the implied volatilities is reasonable, especially between the reconstructed index in 
the simplified model and the index in the original model. 
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Moneyness 

Figure 8: Implied volatility of an individual stock. 



Concerning the market model, by construction, we have the same implied volatility for an individual 
stock as in the simplified model but the implied volatility of the index is far from the simplified one. This 
phenomenon is well known in practice (see P],[5] or [7]) : the implied volatility smile of an index is much 
steeper than the implied volatility smile of an individual stock. The market model of constantly correlated 
local volatility dynamics for the stocks is unable to retrieve the shape of the index smile. A more sophisticated 
dependence structure between stocks is needed. Local correlation models provide an extension of the market 
model in this direction : the correlation at time t between the Brownian motions driving the local volatility 
dynamics of the stocks is a function p(t,I t ) of the index level. But the way this function p influences the 
index smile is not clear at all. Somehow, our models provide another parametrization of the correlation 
structure in which, the function a, that replaces the function p, can be interpreted as the local volatility of 
the index. Yet, the individual stocks can still be properly calibrated. 



Application: Pricing of a worst-of option 

Apart from handling both the index and its composing stocks, our models are also relevant for the widespread 
financial products that are sensitive to correlation in the equity world, such as rainbow options. 

One example of such products is the worst-of performance option whose payout is referenced to the worst 
performer in a basket of shares. For a basket of M shares, the payoff of a call with strike K and maturity 
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Figure 9: Implied volatility of the index. 



T writes ^ min — K j . Our objective is to compare the prices obtained by our model to the prices 

obtained by the market model of constantly correlated stocks. The parameters of the numerical experiment 
are the same as previously and we set the correlation coefficient p such that all the models exhibit the same 
ATM implied volatility for the index. 

The result, as can be seen in Figure flOl is that our prices are always lower than the market model price, 
especially in the money. Hence, a model with a constant correlation coefficient, calibrated in order to fit the 
at the money prices of options written on the index, will always overestimate the risks of worst-of options. 
The reason is that the correlation level needed to fit the at the money prices is very high. Note that the 
prices obtained with the original model and the simplified one are barely distinguishable from each other. 
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5 Conclusion 

In this paper, we have introduced a new model for describing the joint evolution of an index and its 
composing stocks. The idea behind our view is that an index is not only a weighted sum of stocks but 
can also be seen as a market factor that influences their dynamics. In order to have a more tractable 
model, we have studied the limit when the number of underlying stocks goes to infinity and we have shown 
that our model reduces to a local volatility model for the index and to a stochastic volatility model with 
volatility driven by the index for each individual stock. We have discussed calibration issues and proposed a 
simulation-based technique for the calibration of the stock dynamics, which permits us to fit both index and 
stocks smiles. The numerical results obtained on real data for the Eurostoxx index are very encouraging, 
especially for accelerated techniques. We have also compared our models (before and after passing to the 
limit) to a standard market model consisting of local volatility models for the stocks which are constantly 
correlated and we have seen that they lead to a steeper index smile. Finally, when considering the pricing 
of worst-of performance options, which are sensitive to the dependence structure between stocks, we have 
found that our prices are more aggressive than the prices obtained by the standard market model. 

To sum up, we list some properties of our models depending on the options one wishes to handle in the 
Table below 
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Purpose 


Simplified model 


Original model 


Options written on 
-few (J << M) stocks 
-the index. 


-Simulation of a (J + l)-dimensional SDE : 
(I,S\...,S J ). 
-Exact calibration of (S J )i<j<j and I possible. 


-Simulation of an M-dimensional SDE : 
(S 1 ' M ,...,S M ' M ). 

-Exact calibration of (S j ' M )i<j<j possible 
but requires calibration of all the stocks. 
-Approximate calibration of I M . 


Options written on 
-all the stocks 
-the index. 


-Simulation of an (M + l)-dimensional SDE : 
(I,S\...,S M ). 
-Exact calibration of all the stocks possible. 
-Index value : I t = w jS J t . 
-Approximate calibration of I . 


-Simulation of an M-dimensional SDE : 
{S 1 ' M ,...,S M ' M ). 

-Exact calibration of all the stocks possible. 
-Approximate calibration of I M . 



Table 6: Which model to use and when. 
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Appendix 



In order to prove the Theorems 1 and 2, we need the following technical estimation 



Lemma 4 — Under assumption (H^, for all p > 1, one has 



■J- {1 \/|, sup E(\S{ M r)<C p 

0<t<T V 



0<t<T 

w/iere C p = ^ax ^ |S# M | 2p exp ^2r + (2p - l)(max/3 2 + l)K^j pT 
Proof : By Ito's lemma one has 

|^M |2p = |^,M |2p+ /' t |^|2 P((2p)(r _^ )+p(2p _ 1)( ^ a 2 (S;/ M )+ ^ (s55 p/ )))ds 

JO 

+ f{2p)\S^ M \^{P 3 cj{s 1 1?)dB. + Vj (s, Si> M )dWi) 



(19) 



In order to get rid of the stochastic integral, we use a localization technique : let v n be the stopping time 
defined for each n e N by v n := inf{i > 0; \S{' M \ > n}. Then, using (7l[2]), one has 



E 



(l«J 2p ) = I^' M | 2P + E ^ n \Si' M \^{(2p){r + P (2p- l)(/3f<7 2 ( S ,7f ) + r,](s,Si> M ))ds^ 



< 



|S£ A '| 2 " + ((2p)(r - 6 3 )t {r _ Sj >o } + P (2p-l)(fi + l)K 2 b ) J E (|5^5j 2p ) ds 
So, by Gronwall's lemma and the fact that the dividends are nonnegative, 

W<T,E(|^J 2P ) < |^ M | 2p exp((2rp + p(2p~l)(/3 2 + l)At)T) (20) 
Finally, Fatou's lemma permits us to conclude : 

sup E (\S{' M \ 2 ") < \S^ M \ 2 Pe X p((2rp + p(2p- l)(/3 2 + 1)K 2 ) T) . (21) 

□ 



Proof of Theorem 1 : Using the SDEs © and ([5|), one has 



\I t M ~h\ 2p 



J* (If - 7 S ) - jf * f £ 6jWjS^'^ 5l)j ds 

r t ( M \ M t 

/ E PjWjSi' M a(s, 7f ) - /3/^( s , /.) Us s + E Wj / (*, sr w j,/uy |2/ ' 

Jo V=i / j=x Jo 

2PfP ~ 1 Jo^' ~ Is)2PdS + ^ Jo ^^' M ~ Ws j ds 

/•* / M \ 2 P I M rt 

/ V pjWjSi'Vais, 7f ) - /3/ s a( s , 7 a ) d£ fl + V ^ / S^fo Sj^dW- 

Jo J 'i=i ^ 



2p 



2.3 



Hence, using the Burkholder-Davis-Gundy inequality (see Karatzas and Shreve [T^] p. 166), there exists a 
universal positive constant K p such that 

E I sup \I t M - I t \ 2p ] < 4 2p -\a M + b M + cm + d M ) 

\0<t<T J 



where 

JO 

I I M \ 



• a M 



• b M = T 2p ' 



J2S 3 w 3 Sl' M -5I S ] |,/.s 



= ifpTf- 1 jT E f ^ pjWsSWvia, If) - f3I s a{s, I s ) j 1 
d M =K p T p - 1 j\ fast'' 



\r - 

The term aM is the easiest one to handle : 



«m < r 2 ^ T 2p_1 / T E f sup |/f - 4| 2 p) ds. ( 22 ) 

JO \0<u<s J 

Next, using assumption (7l(2) for the first inequality, Holder's inequality for the second and lemma [5] for 
the third, one gets 



-T M M / p \ 

j ° ji=i j p =i \fe=i / 

< KpK* p T p ~ i I E-E(IIO E II(^' M ) 2 ) ds 

Jo h=i j P =i k=i \fe=i / 

T M M p 1 

< K P K% P T P ~ 1 / E-En<( E ((^' M ) 2p )) ? d S 

** o ■ i - i i_ i 



(23) 



ii=i i P =ifc=i 
M 



< K p K 2 b p T p C p \J2 



The same arguments enable us to control the term 6m : 
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rp2p- 



1 [ T (^28 jWj Si> M -5l)j j ds 



< (2T) 2p ^ 1 



ds 




2p\ 


-M 

') 


) 


2p\ 






ds 



E((SI^-SI s ) 2p )d s 



(m \ 2p T 

< 2 2 P- 1 T 2p C p \J2w J \5 J ~8\J +(2T) 2p - 1 6 2p J E 



sup \lf - I u \ 2p ) ds. 

0<u<s 



For the remaining term cm, we will also need the Lipschitz assumption (7l[3]) 



(24) 



pTP -i ^ T E ^ jWj Si' M a(a, If) - pl s a(s, I s ) j j ds 
< 2 2 *- 1 K p TPKl P C p (j^iVjlfo - + 2 2p - 1 K p T p ~ 1 ((3K„) 2p £ E ^sup |/f - J„| 2 ^ 



< 2 2p ~ 1 if p r p - 1 



- E ((/3/fa( S , /f ) - /?J s cx( S , / s )) 2 *>) 



So, combining the inequalities 



e( sup |/ t M -/ t | 2p ) < Co 

\0<t<T J 



and ([25|) , one obtains 



/ / .1/ \" / i/ 

2 



v 



VJ=1 



M 



2/A 



+d / e( sup \I^-I u \ 2 )ds 

JO \Q<u<s J 



with C = 8 2p - 1 TP(TP + K p K 2p )C p and d = 4 2 P-H 22p ~ 1 ^ P rP ~ 1 (/^ ( x) 2p + (ZT) 2 *- 1 ^ + r 2 *> T 2p_1 ). 
Finally, by means of GronwalPs lemma, we conclude that 



M \ I M 

..2 



2p / \ 2p\ 

M ' 



where 
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Proof of Theorem 2 : The proof is similar to the previous one 
\S j t ' M - S{\ 2p < S 2 ^ 1 [ (r - SjfPfP- 1 f (Si' M - Si) 2p ds + 



(si> M Vj ( S ,si- M ) - si Vj ( 3 ,si))dwi 



2p 



+0- 



2/j 



j\si' M o-(s,lM)-Sio-(s,I s ))dB s ^ 



hence, using the Burkholder-Davis-Gundy inequality, there exists a constant K p such that 

E ( sup \St M - Si\A < S 2 ^ 1 f (r - ^) 2p T 2p ^ f E ( sup \S^ M - Sif) ds 

\0<t<T J V Jo \0<«<s / 



+K p T p ~ l / E ((S^- (s, 5i' M ) - 5^ ( S , Sl)f p ) ds 

T \ 

+K p T p - 1 ( 3f£ E ((Si> M a(s, if ) - SMs,I s )) 2p ) d S j 

Using assumption (7tH]), one gets 

/ T E ((Si' M Vj (s, Sf M ) ~ Sivj(s, Sl)) 2p ) ds < K 2p ( T E f sup \S^ M - Si\ 2p ) ds. 

JO JO \0<u<s J 

Finally, by means of lemma [4] and assumptions (7l[2]) and (7-t[3j) , 
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We deduce using Gronwall's lemma 
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where 
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We conclude by Theorem [1] and the sublinearity of the square root function on R + . 
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We now turn to the L 2p -distance between I M and I 

M 



A I 



.M 



i 7 m_ t - |2p = ^ WjS {> M - e^; 

(A/ 
E^i 5 *' 

M M 2p 

^ E- E IK .v 

3i=l j 2p = lk=l 



■2p 



2G 



So, using Holder inequality, one has 
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